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1. Introduction and results 

A Poisson process A'^ = {Nt)t>o with constant intensity A may be regarded as the counting process 
associated with the partial sums (5'n)neNo of ^ sequence {Xk)ke\ti of independent random variables, i.e. 

n 

Nf = sup{n G iNio : Sn < t} with Sq = 0, Sn = for n> 1. 

k=l 

Here each interarrival time has an exponential distribution with parameter A, 

P{Xk <x) = 1 - e"^"^ for all x > 0, 

which we abbreviate to ~ Exp(A). In renewal-theoretic terms, for which see e.g. @], i^t + l)t>o is 
the renewal process associated with the lifetimes (Xfc)fcgiig. 
With 1a the indicator function of the set A, let (/n)neN be defined by 

n-l 



Ii = 1, ^ = n l{^n>x4 for n > 2, 



k=l 

so that In indicates whether Xn is a record within the sequence Xi, X2, X^, ... or not. We are interested 
in the various aspects of the record structure of N, in particular in the two processes C = {Ct)t>o and 
W = {Wt)t>o given by 

Ct ■■= #{l<n<Nt: In = 1}, 
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and 



Wt := [ In,+i ds = Leb({0 < s < t : In,+i = 1}) 
Jo 



for all t > 0; we put = if A'^t = 0. In words: Cj is the number of (completed) record inter-arrival times 
observed up to time t, and Wt is the total amount of time up to time t that the Poisson process spends 
in records. All these variables are defined on some background probability space {^l,A,P). Figure [T] below 
shows these processes on a finite time interval for a particular a; G fi. 
Our first result is a strong law of large numbers for the individual random variables in these processes. 

Theorem 1.1 : With probability 1, and as t ^ oo, 

nm = 1 ana lim = — . 

i^oologf t^oo (logt)^ 2A 

The second part can be seen as an instance of the length biasing phenomenon: if we select an index 
i uniformly at random from {1, . . . ,n}, then the probability that Xi is a record is about (logn)/n if n 
is large. However, if we select t uniformly at random from the time interval [0,n/A] (which contains n 
renewals on average), then the probability that we hit a record renewal is about (logn)^/(2n) for n large. 

A standard question emerging from Theorem 11.11 is whether we can re-scale the respective differences 
such that we obtain non-trivial limit distributions as f ^ oo. Our method for dealing with this question 
gives a considerably stronger result: we show that the whole processes converge. To make this precise we 
use the framework given in the classic 0]. In particular, we regard the processes as random elements of 
the space -D[0, 1] of cadlag functions on the unit interval, endowed with the Skorohod topology and the 

associated Borel structure; denotes convergence in distribution. We write Y(t) instead of Ij if this 
is typographically more appropriate. Then, in order to obtain a converging sequence {Yn)n£i^ of processes 
Yn = {Yn{t))o<t<i from a given process Y = {Y{t))t>o, we put Yn{t) = an{Y{(j){n,t)) - f{n,t)), where, for 
each n € iNi, a,i G R, f{n, •) G D[0, 1], and (pn '■ [0, 1] [0, r„] is a time transformation, by which we mean a 
strictly increasing and continuous function that maps the unit interval onto the initial time segment from 
to r„, where r„ — > cxd as n — > oo. From Theorem 1 1.1 1 it should be clear that the usual scaling (j){n, t) = nt 
would lead to processes with constant paths in some cases; we will instead use (p{n,t) = e"* — 1, so that 
r„ = e" — 1. In particular, we define the re-scaled processes C„ and Wn by 

C^{t) := n-i/2(C7(e"*-l)-nt) 



and 

-^"W -;^(«'(="'-i)-^). 

with < t < 1 throughout. 

In our functional central limit theorem we consider the two processes simultaneously, which means that 
we regard ((?„, Wn) as a random element of the product space D[0, 1] x D[0, 1], endowed with the product 
structures (as is well known, this works for the measurability as well as for the topology aspects because 
of the separability of the Skorohod topology; cf. [1], p. 225). 

Theorem 1.2 : As n ^ oo, 

{Cn,Wn) ^ Z:=(Z^,Z^), 
where the limit process [Z^" , Zj^)o<t<i can be constructed from a standard Brownian motion B = (i?f)o<i<i 
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Figure 1. (a) Poisson process, (b) record indicators, (c) number of records, (d) cumulative time in records 



via 

Zf = -Bt and = [ B^ds- tBt, for < t < I. 

Jo 

Hence, asymptotically and at the level of detail considered here, the second is a simple explicit function 
of the first. 

As will become clear in the proofs, we can replace the discrete parameter n in the definition of the 
processes by a continuous parameter r, say, and the convergence holds with r ^ oo in the sense that we 
have convergence for each sequence (T„)„gy with — > oo as n ^ oo. In particular, we could instead use 
the scaling (j){n,t) = (1 + n)* — 1, which transforms the unit interval into [0,n]. 

We consider the VF-part in more detail. A straightforward calculation (to be given at the end of ^ 
shows that is a centred Gaussian process with covariance function cov{Z^ , Z^^) = min{s^, f^}/3, 
hence Z^ can be seen as a Brownian motion run by a deterministic non-linear clock: the limit process 
Z^ is equal in distribution to [B{t^) / ^y3) ^^^^-^ where B{-) is another standard Brownian motion. 
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Corollary 1.3: As n ^ oo, 



(logn)3/2V ^ ' 2A /o<t<i 

where B = {Bt)o<t<i denotes standard Brownian motion. 

In particular, with the time transformations (j)n(t) = (n+ 1)*^^^ — 1 the initial segments of W have station- 
ary and independent increments in the limit. Of course, Theorem 11.21 also supplies the limit distribution 
for the individual random variables. 



Corollary 1.4: As t ^ oo, 



(logt)3/n ' 2A 

where W' has a normal distribution with mean and variance 1/(3A^). 

Apart from asymptotic normality for the individual variables Wt the theorem can also be used to show 
the asymptotic behaviour of functionals of the process W that depend on the whole history of W up to 
time t, such as in the following corollary. 

Corollary 1.5: As t ^ oo, and with B as in Corollarv \1.3l 



XV3 
[logty/'^ o<s<t 



(log(l + s)) 



2A 



— > sup \Bs\. 

0<s<l 



The limit distribution in Corollary 11.51 is known explicitly; see e.g. 0], p. 80. 

Proofs are given in the next section. In ^we discuss these results and relate them to the literature. We 
also point out some variants and extensions. 



2. Proofs 

We require some more notation. Let 

Li = 1, Ln+i = infj/c > Ln : 1^ = 1} for n > 1 

be the sequence of record times. Then {Rk)k<^ts^ with Ri = Xl, is the subsequence of record values among 
the original sequence {Xi:)^^^ of lifetimes. Further, let 

n n 

An = ^ h, Tn = ^ Rk 

k=l k=l 

be the number of records among the first n lifetimes, and the sum of the first n record values, respectively. 
Then 

is the sum of the record values among the first n lifetimes. We will also need to keep in mind that 
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Proof of Theorem 11.11 : Clearly, 

Zn, <Wt< Zn,+i for all t > 0. (1) 

We need the following two well-known facts from the theory of records of i.i.d. (independent, identically 
distributed) sequences; see e.g. [l] and [4']. First, it follows from the lack of memory property of exponential 
distributions that the sequence (Afc)jtgii^ of differences of successive record values, 

Ai = Ri (= Xi), A„+i = Rn+i - Rn for n > 1, 

is again a sequence of independent random variables with distribution Exp(A). Secondly, it follows from 
the permutation invariance of the joint distribution of the X-values and the continuity of Exp(A) that the 
record indicators In, for n S INJ, are independent and that P{In = 1) = ^/n for all n G INJ. The second fact 
implies that An/\ogn — > 1 a.s. (almost surely) as n — > oo; see Lemma l2 . II b elow . From the first fact and 
the strong law of large numbers we obtain Rn/n ^ 1/A a.s., which together with the elementary fact that 

1 " a 
lim an = a =^ lim } kau = - 

n— >oo n— >oo n-^ ^— ' Z 

k=l 

for real sequences (a„) implies that Tnjv? 1/(2A) a.s. as n — > oo. With probability 1, An T oo, so we 
can pass to the subsequence and obtain 

Zn Ta„ f An I 

a.s. as n — > oo. [Z) 



(logn)2 A^ \\ogn) 2A 

For the second time transformation, from n to Nt, we use that Nt/t A almost surely as t ^ oo; see e.g. 
p. 107 in Q). That implies 

logiV, log(jV,/t) 

— = 1 H ; 1 a.s. as t ^ oo. (3) 

logt logt ^ ^ 

Since Ct = Ai\f^ we gain the first conclusion of the Theorem by replacing n by Nt in the conclusion of 
Lemma [2Tt and applying ([3]). 
Finally, from ([T|) we obtain 

Zn, flogNt^-^ ^ Wt ^ Zn,+i flogiNt + 1)^^ f log Nt^"^ 



(log iVt)n log t ; - (logt)2 - (iog(Ar, + i))H log Nt J \logt 

The result for Wt now follows on using ([2]) and ([3D, together with the elementary properties of almost-sure 
convergence. □ 

Lemma 2.1: j4„/logn 1 a.s. as n ^ oo. 

Proof: Let Bn = {In - n"^)/log(n + 1), then the B n are independent, of zero mean, and E{Bl) = 
{n~^ — n~2)/(iog(n -|- 1))^; consequently Yl'^=i -^i-^n) < It follows, for instance by Theorem 
12.2(a), that X^^i converges almost surely. By Kronecker's Lemma, Yl^=ii^k — k~^) / log{n -|- 1) ^ 
a.s. The claimed result follows. □ 

For the proof of Theorem 11.21 we first observe that by a straightforward re-scaling argument, we may 
assume A = 1. We keep that assumption in force up to the end of Section [2j To start the proof, note that 

Wt = Ta,^ + {t-SNjiN.+l. (4) 
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The distribution of the 'current age' t — S^t is known exphcitly for Poisson processes (it is a truncated 
exponential) and generally converges as t — > oo as long as the lifetime distribution has a finite mean. In 
particular, we have {t — SNt)lNt+i = Op(l), and the process version of this bound will allow us to ignore 
the second term on the right hand side in It is straightforward to obtain asymptotic normality for 
(Tn}n£ts^, but the time change (from n to A^J requires work. Anscombe's theorem (see e.g. Theorem 7.3.2 
in [sl]) requires zero-mean summands and then leads to asymptotic normality for {Ta,^ — A^/2)/(log n)'^/^ 
only. This cannot be combined with a similar statement for An unless we have the joint distribution. Hence 
the proof below provides the joint distribution together with the time change(s). It proceeds by a sequence 
of lemmas. 
We define 

[x\ = max{n £ Z : n < x}, [x] = [x\ + 1, [x] = min{n G Z : n > x}. 
The first two of these are right-continuous. 

Lemma 2.2: The sequence {Rn) of record values satisfies afunctional central limit theorem: 

Rn := {n-^/\R[nt] - ni))o<t<i B = {Bt)o<t<i, (5) 
where B denotes Brownian motion. 

Proof : We work with the results in Section 17 of [1]. To obtain dS}, it is enough to invoke Donsker's 
theorem together with the result on the record-value differences that was already used in the proof of 
Theorem 11.11 □ 



Lemma 2.3: We have 



n 



-'/^n-%^t^- I R^nu]du]) ^0, (6) 

0<t<l 



and 

{n-'/\\ogL[nt]-Rint]))o^t^, ^ 0. (7) 

Proof: Note that the union of the disjoint intervals ((i — l)/n,i/n] for i = 1, 2, . . . , [nt] is the interval 
(0, i] plus the interval (t, Adding the integrals on the intervals ((i — l)/n,i/n], and subtracting 

that on {t, [nt]/n], we obtain 

nt] 

nt] 



where 9 = 9nt lies between —1 and 0. Thus 



n ^^"^ sup 

0<t<l 



n ^T^nt] - / R[nu] du 



t 

< n"^/^ sup \Rint]\ 

0<t<l 

< sup \Rnit)\ +n~^/^ (8) 

0<t<l 
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By ([5]) the supremum converges in distribution to supo<«x \Bt\- So the right-hand side of dS]) converges 
to in probabihty, which proves 

For ([7|) we start from the representation for the record time process given in 0, Theorem 5.8. With 
(^)ieN another sequence of independent standard exponentials, also independent of (i?„)„eiH, we have 



EN 



where 



n-l 



1=1 



log(l 



-Ri 



and '=' denotes equality in distribution. 

Set Ml := min{n : i?„ > log 2}. For i > Mi we have e~^' < 1/2, so the fact that x < — log(l — x) < 2x 
for < a; < 1/2 gives e~^' < - log(l - e"^') < 2e~^\ So for n > Mi, 



n-l 

i=Mi 



log(l 



n-l 

4f! + E ^^^^^^ 

i=Mi 



Ri 



where 1/2 < (pi <1. 

Now let M2 be the least m > Mi such that T^e^* > 1 for alH > m. We need the fact that V^e^* tends 
to infinity a.s., which is left as an exercise for the reader. It implies that M2 is a proper (finite) random 
variable. For i > M2 we have Vi(f>ie^^ > 1/2, so the effect of the rounding up given by the [ and ] brackets 
can be considered as multiplication by ipi where 1 < < 2. Thus for n > M2, 

n-l 



i=M2 



By setting (pi = ipi := 1 for i < M2 and defining Ln^ = 1 + J2i=i Vi4'i'^i^^' foi' all n > 1, we obtain 

lW=A + l(,2) forn>M2, 



where A := L 
Next, 



(1) 

M2 



aRi 



logL« =logLl2) + log 1 + 



A 



L 



(2) 



for n > M2. 



For the last term here we use the fact that Ln > L^^l for n > M2, which implies that | log(l + A/Ln 



(2) 



-(2)^ 



('2) z 

increases when we replace Ln in it by L^^^ (consider the cases ^ > and j4 < separately). For the 

other term on the right-hand side we use the fact that 1/2 < (piipi < 2 to deduce that log Ln differs from 
log(l + YJlzl Vie^^) by at most ±log2. Therefore 



'(2) 



n-l 



(logL„ - Rn)nm =Un + log(e-^" + ^ yie-(«"-«') 



i=l 



Charles M. Goldie and Rudolf Griihel 



where 



|en| <log2 + 



log 1 + 



A 



L 



for n > M2. 



It is clear that n sup„gy \en\ > as n ^ 00. Formula ([7]) follows by virtue of Lemma 12.41 below. □ 
Lemma 2.4: With the Vi as in the previous proof, 



n ^^"^ max 

k=l,...,n 



fe-1 



i=l 



as n ^ 00. 

Proof : Choose e > 0, then 



/ k-l \ fe-l 

P log(e-«'' + J2 ^ie-^^---^')) >e^k \ <p{l + Y,V^> 
\ 1=1 / i=i 



E{1 + E-=l V^ 



< 



As the latter is summable it follows that, with probability 1, there are at most finitely many k such that 
log(e^-^*^ + ^^^Z\ T/je~*^^'=~^')) exceeds e-^/A;, and it follows that 



fc-i 



max log(e-«^ + ^l/ie-(^'=-«')) -^0. 



(9) 



On the other hand, 



fc-i 



fc-i 



P log(e-^^ + ^F,e-(^^-^-) 



i=l 



= P{Vu-i < e-'^'') + P{Ak >eVk + log V^-i > 0) 

= 1 - exp(-e-^v^^) + r e-(^^+'°s^')e-^dx. (10) 



In the last right-hand side the term 1 — exp(— e '^^^) is (1 + o(l))e "^^^ as A; — > 00. And the integral is the 
sum of 



-ey/k 



dx 



X 



{Ei{-) is the exponential integral), and 



-e^k 



dx 



.^fc X 



dx 



e^k X 



Thus the left-hand side of (jlOp is bounded by a constant multiple of ^/k e and so is summable. Arguing 



Cumulative record times in a Poisson process 

as for Q, it follows that 

fc-i 



n"^/^ max -logle'^" Vie'^^'-'^^A ^ 

k=l,...,n V ' / 



and with ([9]) this gives the result. □ 
Define 

fn:= (n-3/2(r[„,]-(nt)V2))o<,<,, 

the (standardised) record sum process, and 

L„ := (n"^/2(logL[„i] -nt))^^^^^, 

the (standardised) record times process. 
Lemma 2.5: As n — > oo, 

where we define ^ : Z)[0, 1] ^ L»[0, 1] x L>[0, 1] 6y ^'(/) := (^'i(/), /), with *i : L»[0, 1] ^ ^[0, 1] in turn 
defined by ^i{f){t) = j'^ f{s) ds. 



Proof: Clearly, is continuous (and linear), so that ^ is continuous too. By the continuous mapping 

theorem, applied to the conclusion of Lemma 12.21 it follows that '^{Rn) — -> '^{B). But Lemma 12.31 savs 
that {Tn,Ln) = ^(-Rn) + op{l), SO our claimed result follows. □ 

Now we switch to the 'stochastic clock(s)' given by 

$„(t) := -^jv(e"*-i), for < i < 1, 

and deduce joint convergence of the time-changed processes T„ o and L„ o <!)„. 
Lemma 2.6: As n oo, 

Proof : From Lemma 12.11 it follows that 

max |j4fc — log(l + /i;)| = op(n). (11) 

l<fc<e" 
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Choose e > 0. Then 



P sup |^7v(e--i) - log(l + A^(e"* - 1))| > en 

0<t<l 



< P{N{e'' - 1) > e"+i) + P 7V(e" - 1) < e"+\ sup - log(l + iV(e"* - 1))| > en 



iAr(e"t_i) 

0<t<l 

< P(7V(e" - 1) > e"+^) + P ( 7V(e" - 1) < e"+\ max | - log(l + A;)| > en 

l<A;<e"+i 

< P(7V(e" - 1) > e"+i) + P ( max - log(l + > -(n + 1) ) , 

^l<fe<e"+i 2 

since en > e(n + l)/2. On the right-hand side, (jlip makes the second term tend to 0, while the first term 
tends to because Nt/t — > 1 a.s. We have shown that 

sup - n-^ log(l + iV(e"* - 1))| = op(l), 

0<i<l 

and it follows immediately that 

sup -t| = op(l). (12) 

0<t<l 

Using this, we may proceed as in the proof of Theorem 17.3 in \§\ to conclude the present proof. □ 

Inserting A]\r^ for n in L„ is like time inversion, so we expect that log L^^^(i~^ differs from nt by a negligible 
amount. That is indeed the conclusion of Lemma |2 . 1 2 1 b elow . and the intervening results are steps on the 
road to that end. We allow ourselves to write At for non-integer t, to be interpreted as ^i^^j , and similarly 
for other random sequences. 



Lemma 2.7: As n ^ oo, 



sup n-^/^logL{Ae^t + 1) - logL(^e-)) = op{l). 

0<t<l 



Proof : For m > k, 



so 



k k+1 m-1 k 

P(L„+i > m\Ln = k) = -r—T—^ = — , 

k + 1 k + 2 m m 



k 1 

P{Ln+l > xLn\Ln = k) = -—— < 



which we use in the form 



\xk\ X — 1 ' 



P\^>x\<^ {x>l). 

Ln / X - 1 
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Choose e > 0, then 

P ( sup ^^fr; > e^^A < P{A,. >2n)+p( max > e^^A 

< P{A,. >2n) + yp(^> e'^A 



2n 

< P{Ae^ > 2n) + 



eV" _ 1 ' 

and this tends to as n ^ oo because Ae^/n — > 1 almost surely. This establishes Lemma 12.71 □ 
Lemma 2.8: As n oo, 

(n-i/2(logL(^e-)-ni)) =op{l). 
\ ^ Vo<t<i 

Proof: For < i < L^+i, Ai = k, so = Lj^. Therefore L^^ <t< L^t+i for all t. We use this in the 
form 

< logt - log La, < logLA,+i - log La,. 

Thus 

sup n"^/^|logL(Ae«t) -nt| < sup n"^/^ (log L(Ae"' + 1) - log L(^e-.t)) , 

0<t<l 0<t<l 

and the result follows by Lemma 12.71 □ 
Lemma 2.9: As n ^ oo, 



n 



^1/2 



(logL(^e-_i + l)-nt)) =op(l). 



0<t<l 



Proof : We write 

"^/^(logL(Ae,.t_i + 1) - nt) <( sup + sup ) n"^/^(logL(^e"'-i + 1) - nt) 



sup 

0<t<l 



n 



0<i<n-ilog2 n-ilog2<t<l 



The first supremum on the right equals 



sup n ^/^l logLi - nt| = n log 2 = o(l). 

0<t<n-i log 2 

The second is bounded by 

sup n"^/2|logL(Ae.t_i + 1) - log(e"* - 1)1 + sup log(e"* - 1) - nt| 

n-ilog2<t<l n-ilog2<t<l 



< sup n ^/^|logL(Ae"' + 1) - nt| + sup n ^/^|log(l 

0<t<l n-ilog2<t<l 



< sup n"^/^|logL(Ae"f + 1) - logL(^e,.t)| + sup n"^/^|logL(^e«t) -nil +o(l). 

0<t<l 0<i<l 

The two terms on the final right-hand side are op(l) by Lemmas 12.71 and 12.81 respectively. □ 
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Lemma 2.10: As n ^ oo, 

sup n-i/2| iog(Ar(e'^* _ i) + i) _ nt\ = o(l) a.s.. (13) 

0<t<l 

Proof: From Nt ^ t a.s., i.e. Nt/t 1, it follows that N{e* — 1) + 1 ~ e* a.s. Therefore 

M := sup I log(iV(e* - 1) + 1) - t| 

0<t<oo 

is a proper (finite) random variable. The left-hand side of is bounded by n~^/^M, so converges to 
a.s. as n ^ oo. □ 



Lemma 2.11: As n ^ oo, 



n~^/^{logL{ANie-t-i) + 1) - nt))^^^^^ = op{l)- 

Proof: First, because of Lemma I2.1UI it suffices to show 

n -V2(iogL(^^(^„,_,) + 1) - log(iV(e-* - 1) + 1)))^^^^^ = op{l). 

Choose e > 0. We have 

pf sup n-i/2(iogL(^^(e.t_i) + l)-log(iV(e"*-l) + l)) > 
Vo<t<i / 

< P(iV(e" - 1) + 1 > e^") 

+ P fiV(e" - 1) + 1 < e^", sup n-i/2(iogL(^;v(e--i) + 1) - log(iV(e"* - 1) + 1)) > e] . 
V o<t<i / 

Set s{t) := n^^ log(A^(e"'* — 1) + 1), then as t varies over its range < t < 1, s varies over the range 
< s < n^"'^ log(A^(e"' — 1) + 1), which is contained in < s < 2 because we are within the event 
N{e^ — 1) + 1 < e^"'. Thus the above is bounded by 

P(iV(e" - 1) + 1 > e^") + P ( sup n"^/2|logL(Ae.= _i + 1) -ns| > e ) 

Vo<s<2 / 

= P(A^(e" - 1) + 1 > e^") + P ( sup 2^^^m-^/'^\logL{Ae^t^i + 1) - mt\ > e) , 

where we set t := s/2 and m := 2n. In the final right-hand side the first term tends to because Nt/t — > 1 
a.s., and the second does so too because of Lemma \T9\ □ 



Lemma 2.12: As n ^ oo, 

^nO^n = -C„ + Op(l). (14) 

Proof : Observe (recalling that [k] = k + 1 for integer k) that 

Ln o^n= (n"l/2(logL(A^(e„t_l) + 1) - ^TVCe-'-l) 

= (n^^/'^{logL{AN(^e^.t_i^ + 1) - nO)^^^^^ + n^^\t - <^n{t))o<t<i- 



0<t<l 
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By Lemma Pi. Ill the first term on the right is op(l) as n ^ cxd. Since Cn{t) = n^^'^(^^n{t) — t), we conclude 
{[HD. □ 



Proof of Theorem 11.21 : From Lemma 12.121 we conclude that 

(r„ o Cn) = {fn O -L„ O $„) + Op(l), 

and then from Lemma 12.61 that 



{fnO<^^,Cn) ^(.^l(.B),-B). (15) 



Now 



^ ^ ' J ^<t<\ \ ^ ^ / z /0<t<l 

= r„ o - [\n^l'^ {t - <Pn{t)) {t + ^n{t))) 

\^ / 0<t<l 

= rnO$„- (nl/2(t-$„(t))t) +Op(l), 

V ^ /o<t<i 

with the last step by ()12p . From (llSp . as addition is continuous in D[0, 1] x D[0, 1] at elements (x,y) with 
no discontinuities in common (^13], §4), it follows that 

(^Cn, (n-3/2(r„^^(,) - i(nt)2)^^^^^J ^ (-5,^1(5) - K)o<i<i) • (16) 

Now from (jlj), writing t := e"* — 1, 

W{t)=TA,^,^+{i-SN,)lN,+l 

= Tn<S,„{t) + {i- SNf)lNt+l, 

SO that 

Wn{t) = n-3/2(r„$^(,) _ i(nt)2) + n-3/2(t _ 5^J/^^.^,. 

We will thus have 

Wn = n-^/^(T^^^^,)-hntf) + op{l) (17) 

provided we can show that 

-3/2 ,. P 



n sup {t-SNt)lNt+i — ^ 0. 

0<i<e"-l 



We prove the equivalent assertion 



(logn)-3/2 sup (t - 5jvJ/jv.+i 0. (18) 

Q<t<n 
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Choose e > 0, then 



P((logn) sup {t - Sn,)In,+i >e)< PiNn > 2n) + P((logn)-3/2 max{Xi, . . . , X2n+i} > e) 

^ 0<t<n ' 



(logn)-3/2\2n+l 



= P{Nn > 2n) + 1 - (1 - e-^('°s")'''') 
2n 

= o(l) + (l + o(l))-^(^^^0. 

This proves ([TH]) . 

Combinmg (I16p and ()17|) . we have thus proved that 

{Cr,,Wn)^{Z^,Z'^) with Z^ = -Bt and = ^-1(5)* - tS^. 

It remains to find the covariance function for the hmit process . Let Yj = ^'i(i?)t = f^Bgds. For 
s < t we obtain 

EYgYt = E I I BuB^dudv= f f mm{u,v}dudv = —s^t s^, 
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which together with the results 



E{BsYt) = st- h^, E{BtY,) = 



similarly derived, leads to 



cov(Zf,Zf) = EYsYt-sEB.Yt-tEBtYs + stEBsBt = 



□ 



3. Complements and remarks 



3.1. 



A weak version of Theore mll.il with convergence in probability rather than almost sure convergence, has 
already been obtained in [iz 



3.2. 

For a lifetime distribution with bounded support and an atom at the supremum of its support the number 
of record renewals is obviously finite with probability 1. For discrete uniform distributions on {1,2,... , d} 
asymptotic normality was obtained in [13] for the total time spent in records as d — > 00. That paper also 
contains some remarks on the interface of renewal theory and the theory of records, manifesting itself 
there in the 'perpetuity' Yl^=i \Yk=i with {Ukjk&Vi a sequence of independent random variables, each 
uniformly distributed on (0, 1). A perpetuity appears in the present paper also, at the end of the proof of 
Lemma 12.31 and in the Lemma following, in the shape of the sequence 

n-l 

1=1 
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This propagates itself by iteration of the relation Yn+i = (Yn + Vn)Un+i, where Yn, Vn and Un+i are 
mutually independent, Vn has a unit exponential distribution and Un+i = e^^"+^ is uniform on (0,1). The 
stationary distribution is that of y = Yl'^=i ^ 11^=1 ^k- Perpetuities of this type are discussed in [8[. The 
results of Lemma 12.31 can be seen as remarkable in that they show that the record value sequence takes 
most of the randomness out of the record times sequence; the randomness that 'remains' is encapsulated 
in the perpetuity. 

3.3. 

We could divide the time line into record and non-record pieces, as we have, not just for Poisson processes 
but for completely general renewal processes. How does the tail behaviour of the lifetime distribution 
affect the relative size of the record pieces? Above we have worked out the answer for exponential lifetime 
distributions, but we expect it to be possible to find analogues of Theorem 11.11 for other distributions. 
The 'exponential quantile function' (i.e. the monotone function that turns a standard exponential into a 
random variable with prescribed distribution) is likely to play a key role. 

3.4. 

Our results can also be seen in the larger context of ordinary sampling vs. renewal sampling, and its 
consequences for the theory of records. In order to explain this we start with a one-dimensional distribution 
function F which, for simplicity, we assume to be continuous. 

By an ordinary sample of size n we mean a family Xi , . . . , X„ of independent random variables with 
distribution function F. Like many of the core results of probability theory the classical theory of records 
refers to this setup. In applications, however, the random variables Xi, X2, . . . , while still being inde- 
pendent and identically distributed, often arrive as the times between successive events in some counting 
process, so that instead of Xi, . . . , X„ we observe the associated counting process N = {Nt)t>o over a 
fixed finite time interval. Obviously, this situation, which is what we mean by renewal sampling, changes 
the stochastic structure of the data. It is trivially true, for example, that observing over a fixed time 
interval [0, t] has a truncation effect as the values of F then matter only on this interval. Also, and more 
interestingly, the celebrated inspection paradox can be seen in this context: it essentially means that the 
distribution C{Xn^^i) of the inter-arrival time straddling t differs from the original distribution of the 
A- variables. 

How does renewal sampling change the stochastic structure of records? One of the felicitous aspects 
of the ordinary sampling case is the fact that for many results the distribution of the A-variables does 
not matter (as long as F is continuous). That may no longer be true if we replace ordinary sampling by 
renewal sampling. Indeed, with 

Cn := #{l<i<n: Xi> A„} 

the rank of the last variable in the sample, it is known that C„ is uniformly distributed on the set 
{1, . . . , n}, so that the distribution of n~^C„ converges to the uniform distribution on the unit interval as 
n 00, whereas it was shown in [§] that this is no longer true if we consider {Nt + l)~^C]\ft+i instead. 
For example, in the case of a unit-rate Poisson process, i.e. with C{Xi) exponential with mean 1, we still 
have convergence in distribution, but the limit distribution now has density x 1— > — logx for < x < 1. 

Of course, not all aspects of the record structure will change when passing from ordinary to renewal 
sampling. Consider for example the maximum 

Mn := max{Ai, . . . 

of the first n values for samples from the exponential distribution with mean 1. It is well known that 
Mn — logn converges in distribution as n ^ 00, and that the limit distribution is the Gumbel distribution, 
and it is straightforward to prove that the same holds for M^^ — log Nt as t 00. Nevertheless even 
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this example can serve as the basis for a seemingly paradoxical consequence of random sub-sampling: the 
current maximum of the completed lifetimes is equal to the current record value, so that M^t = Ran^ 
the notation used in the previous sections. Now we know that {Rn)n£ts4 is the sequence of partial sums 
associated with an i.i.d. sequence of standard exponentials, which implies that {Rn — n)/^/n converges in 
distribution to a standard normal as n ^ oo. Thus, sampling along A^^ removes the need to scale the 
centred variables by a factor that tends to infinity, as for the maximum a simple shift is enough to obtain 
a nontrivial limit distribution (of course, the point is that the shift obtained by inserting Nt for n is now 
random) . 

Finally, we wish to point out that the behaviour of records under renewal sampling has applications 



in the analysis of algorithms. In a discrete variant of renewal sampling, with geometric inter- arrival 
times, appeared in connection with the analysis of an iterative procedure, von Neumann addition. There 
the length of the longest time interval between successive events essentially determines the number of 
iterations needed, which in turn is the main influence for the running time of the algorithm. 
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